Capillary condensation revisited 
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We present new results and insights on liquid-vapor phase transitions ol a fluid with long-range 
dispersion intermolecular forces spatially confined by a substrate that forms a rectangular cavity 
and also exerts a long-range dispersion potential. Our analysis, based on the statistical mechanics 
of liquids, namely density functional theory, allows for a systematic and complete exploration of the 
system's phase space, in particular density configurations and phase diagrams. 
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Capillary wetting phenomena are crucial in a wide va- 
riety of natural processes and technological applications. 
From the wetting properties of plant leaves [1] to build- 
ing of nanoreactors [5] and design of "lab-on-chip" de- 
vices [3]. At the same time, they provide one of the most 
striking manifestations of the attractive intermolecular 
forces governing the behavior of matter. Not surpris- 
ingly, there has been a vigorous interest in such systems 
for several decades generating an abundance of funda- 
mental research. 

Of particular interest is the phase behavior of confined 
fluids. For microscopic systems concepts such as surface 
tensions and contact angles become quite limited. One 
must explicitly take into account the molecular struc- 
ture of the fluid, since most of the exciting effects are 
caused by the interplay of different length scales, such as 
the ranges of intermolecular potentials and the charac- 
teristic dimensions of the confining geometry. These ad- 
ditional parameters act as independent thermodynamic 
fields, and can lead, according to the Gibbs rule, to a rich 
phase behavior. 

One of the most widely studied phase transitions in 
inhomogeneous fluids is capillary condensation (CC) in 
which a confined space becomes filled with condensed liq- 
uid from the vapor. The main theoretical approaches to 
CC are mean field Van der Waals loops and renormaliza- 
tion group theory. While the latter has been successfully 
applied to estimating the extent of fluctuation effects in 
various confined systems [4j[5], far less attention has been 
given to the molecular structure of the adsorbent. Here 
we adopt a fully microscopic approach based on density 
functional (DF) theory leading to a molecular model for 
the fluid which is not restricted to complete wetting as 
previous works [5]- A two-dimensional (2D) square cap- 
illary is used as a model system. We undertake a system- 
atic investigation of CC and obtain detailed information 
on all basic characteristics of the fluid, including density 
profiles and menisci shapes, along with bifurcation dia- 
grams of possible states as well as coexistence curves and 
saddle points of the free energy. 

Consider a long-ranged interacting fluid confined in a 
rectangular domain of width H and length L by a sub- 



strate whose walls exert long-ranged forces on the fluid. 
The density is assumed constant in the direction orthog- 
onal to the H-L plane: p(r) = p(x,y). The system is 
assumed connected to a thermostat fixing the values of 
the chemical potential p and temperature T. Classical 
DF theory allows us to express the free energy of the 
pairwise interacting fluid in an external field formed by 
the substrate walls as a unique functional of the one-body 
density. Specifically, the grand canonical free energy is 
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where the canonical free energy functional F follows from 
perturbation around a hard sphere fluid in powers of the 
attractive potential <^ a ttr up to first order [5J: 

F [p (r)] = fdv (/ id (p) + pvj (p) + P V cxt ) 



+ -J drj dr'p (r) p (r') p attI (|r - r'|) , (2) 

/id being the ideal free energy. For hard spheres of diam- 
eter a the configurational part of the free energy is [7] 

1>(p) = k B T r,( ?~ 2 $ , r, = *o* P fa (3) 
(1-77) 

the long-ranged potential of the fluid-fluid interactions 
with strength parameter e is [8] 
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and of the fluid-wall interactions 

V(x,y) = Vo(x) + V (y) + Vo(L-x) + V (H-y), (5) 
where 

1 I er \ If a 



V (z) = E -- 



6 V z + z 



45 V z + z a 



. (6) 



2 



defines a wall potential with the parameters of strength 
Eq, range 00 and low- 2; cut-off zq. The density configu- 
ration p(x,y) minimizing D, is obtained numerically via 
a 2D extension of the one-dimensional (ID) numerical 
method developed in [5] , which also allows to obtain (via 
arc-length continuation) adsorption isotherms and phase 
diagrams of the fluid, as well as to assess the thermody- 
namic stability of particular configurations. 

The DF ([2| is known to capture the physics of fluids at 
vapor and liquid densities. The use of ab initio substrate 
potentials in place of ^ [TU] , as well as a different version 
of DF including the repulsive effects of molecular packing 
non- locally in p (r) , through an improvement of the term 
tp(p(r)) in |2]) 11], would possibly allow for comparison 
with simulations or experiments. It would not, however, 
add to the understanding of the physical phenomena de- 
scribed here which are adequately captured by (j2j). 

Consider, first, an open capillary pore, i.e. L — > 00. 
When it is formed by two identical plates (of infinite 
area), a simple thermodynamic expression for the free 
energy required to form two liquid-wall interfaces leads 
to the Kelvin shift of the bulk coexistence curve: 
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where a\ v is the liquid-vapor surface tension and O is 
the contact angle. The pressures P vap and Pa q refer to 
the capillary- vapor and -liquid phases, which, due to Q, 
can form inside the pore from fluid whose bulk is vapor. 
There exists an abundant amount of earlier literature on 
the wetting properties of an infinite ID slit pore, where 
the fluid density varies in the direction orthogonal to both 
walls, the key finding being the discontinuous first-order 
transition to CC |12j . 

A semi-infinite slit pore capped from one end with a 
third wall, on the other hand, is a prototypical 2D capil- 
lary system. It benefits from the results for ID slit pore 
in that the state of the fluid in the latter acts as a spec- 
tator phase for the former. The discontinuous CC of a 
slit pore thus becomes a continuous second-order phase 
transition in the capped capillary, analogous to the com- 
plete wetting of a single planar wall, but corresponding in 
the case of the capped capillary to the propagation of the 
liquid meniscus away from the capping wall as the fluid 
pressure approaches that of CC j4] . Can the analogy with 
the ID wetting be extended beyond that observation? 

Our computations indicate the presence of a discontin- 
uous first-order transition corresponding to the formation 
of the liquid meniscus at a finite distance from the cap- 
ping wall, as the fluid's chemical potential (pressure) is 
increased towards CC. Depending on the wetting prop- 
erties of the walls, the transition can take place between 
density profiles of varying topology as shown on Fig. [T] 
For illustration purposes, throughout the paper we have 
rescaled the 2D data in each case between capillary p vap 
and pnq and defined the interface (dashed curves) as a 
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FIG. 1. Coexisting density configurations in a capped capil- 
lary of width 30er. In each case the data is rescaled between 
capillary p V a P and pn q . The dashed line marks the liquid- 
vapor interface, (a) and (b): all walls are the same and hy- 
drophobic; (c) and (d): top and bottom walls are the same 
and wet, hydrophobic left wall; (e) and (f): bottom and left 
walls are the same and hydrophilic, hydrophobic top wall. 



density level set equal to (pn q + p vap ) /2. We have fur- 
ther limited the description of our choice of parameters 
to providing the wetting properties of the walls during a 
planar adsorption at bulk coexistence as wet (O = 0°), 
hydrophilic (0 < 90°) and hydrophobic (0 > 90°). See 
Movies 1 and 2 of the Supplemental Material for a demon- 
stration of continuous and discontinuous phase transi- 
tions in the semi-infinite capped capillary. 

Our first finding is that wetting in the capped capillary, 
including the discontinuous meniscus transitions (e.g., 
cases in Fig. [lja)-(f)), maps on the ID wetting of a single 
planar wall immersed in vapor, if the shift in the chemical 
potential is counted from the ID slit pore CC at the 
same temperature. Figure [2]ja) shows the typical phase 
diagram of a capped capillary in the Ap - T plane, where 
Ap is counted from the bulk liquid-vapor coexistence. 

The CC line Ap c (T) (full line in the figure), which 
completely describes the wetting of a ID slit pore, plays 
in the case of a 2D capped capillary the same role as is 
played by the bulk coexistence line (Ap = 0) in the ID 
planar adsorption on a single wall: it denotes the limiting 
values of Ap, at which the spectator phase is expected 
to transform, necessitating the continuous second-order 
phase transition of the fluid in contact with the sub- 
strate. Analogously to the pre-wetting transitions dur- 
ing the ID adsorption on planar substrates [9], tangent 
to the Ap c (T) we find the A/i m (T)-line of first-order 
transitions corresponding to the discontinuous formation 
of the liquid meniscus in the capped capillary at a finite 
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FIG. 2. (a): Capped capillary phase diagram. Full line: 
A/i c (T) - ID slit pore condensation line. Dashed line tangent 
to Au c (T): A/i m (T) - meniscus transition line, denoting the 
location of first-order transitions corresponding to the forma- 
tion of liquid meniscus. Dashed line between Td ci and Td C2 : 
A/id (T) - droplet transition line of the first-order transitions 
corresponding to the formation of capillary-liquid droplets in 
the corners, (b): isotherm of the route designated in (a), ver- 
tical dotted lines correspond to A/id and A/i m - intersections 
of the route with the lines of first-order transitions, (c)-(e): 
typical coexisting density profiles in a capillary with all walls 
the same and hydrophilic (see text). 

distance from the wall (dashed line). 

Central to the understanding of the system's surface 
phase behavior are the critical temperatures T cw and T mc . 
For temperatures below the capillary wetting tempera- 
ture, T cw , the values of A/i can isothermally cross the 
CC line (and go up to its spinodal, not shown in the fig- 
ure), but the fluid would remain in the capillary- vapor 
phase (capillary-liquid phase being metastable). This 
phenomenon is completely analogous to the partial wet- 
ting of a single planar wall, when at T < T w and A/i = 
the fluid does not adsorb a liquid layer on the wall and re- 
mains in the vapor phase (0 > 0° at liquid-vapor coexis- 
tence). Further on, in the capped capillary we find a dis- 
continuous jump to CC at T = T cw and A/i = A/i c (T cw ): 
the capillary becomes filled with the capillary-liquid, just 
as a single planar wall becomes completely wet by liq- 
uid at T = T w and A/i = 0. Obviously, as H — > oo, 
Tcw ~ > T w of the capping wall, and since the value of 
T cw is determined by the interplay of all relevant ther- 
modynamic fields in the system (strengths and ranges of 
potentials, distance of wall separation, as well as the bulk 
fields T and /i), the dimensionless ratio fee (T cw — T w ) je 
can serve as a measure of the confinement- specific capil- 
lary effects in possible experimental realizations. 

In the range of temperatures between T cw and the 
meniscus critical temperature, T mc , the fluid exhibits a 
discontinuous first-order meniscus transition, analogous 
to pre- wetting in the case of planar adsorption [5] . Above 



T mc , the liquid meniscus is being formed continuously as 
A/i — > A/i c (T) isothermally. 

Figure [2^b) depicts the typical adsorption isotherm 
corresponding to the route designated on the phase dia- 
gram, Fig. [2]ja). The excess grand free energy f2 ex , de- 
fined as the difference between the total grand free en- 
ergies of the fluid inside the capped capillary and the 
ID slit pore of width H at the same values of T and 
/x, possesses the typical hysteresis features of the mean 
field Van-der-Waals loops. The fluid configurations coex- 
isting at A/x = A/i m are capillary- vapor (Fig. [2jc)) and 
capillary- liquid slab with a meniscus (Fig. [2](d)). The in- 
set on Fig.[2jb) shows an additional first-order transition, 
at A/i = A/id, between two metastable configurations of 
capillary-vapor, and capillary-liquid drops formed in the 
corners (Fig. [2je)) , demonstrating the interplay between 
the capillary- and wedge-filling phenomena. The phase 
line of that latter transition, A/id {T) (thin dashed line in 
Fig. [2]ja)), terminates at the critical temperatures Td Cl 
and T dC2 . 

For isothermal routes A/i — > A/i c (T), with T < T^ Cl 
droplets do not form in the capillary corners, the transi- 
tion to the phase with meniscus (e.g., Fig.[2](d)) happens 
from the vapor phase (e.g., Fig. [2^c)). For T > Td C2 the 
formation of droplets in the corners happens continuously 
from the vapor phase as A/i — ► A/i c (T), the meniscus 
transition takes place between the configurations analo- 
gous to those in Fig. [2^e) and Fig. [2^d) . Increasing H will 
lead to higher isolation of wedges in the corners of the 
capillary, shifting the vapor-to-droplet line down in the 
phase diagram, and allowing for a triple-point, where all 
three phases of types shown on Figs. [2jc)-(e) coexist, sug- 
gesting also the possibility for the existence of a critical 
value H = Ho > 1, such that for H > H Q and T < T cw 
the complete saturation of capillary with capillary-liquid 
is achieved via the merging of the two corner droplets, 
without a meniscus forming (see, e.g. Fig.[ljc)). 

Our second finding is related to the second-order phase 
transition corresponding to the critical unfolding of the 
meniscus as A/i — > A/i c (T). The earlier work based on 
the effective Hamiltonian model for the capped capillary 
by Parry et al. [I] offered the critical exponents for the 
divergence of the order parameter, which was the dis- 
tance of the liquid meniscus to the capping wall. On the 
other hand, our fully microscopic approach is based on 
exactly accounting for the intermolecular effects and al- 
lows us to consider the natural physical order parameter 
of the system, which for A/i — > A/i c (T) is the thermody- 
namic variable conjugate to the chemical potential, i.e. 
the adsorption: 

r = j dr(p(x,y)-p 1D (y)), (8) 

where pid (y) is the fluid density inside the ID slit pore 
serving as the spectator phase, the integration is done 
over the entire volume of the infinite capped capillary. 
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FIG. 3. (a): Typical adsorption isotherms of a capped cap- 
illary. In the presented case H = 30, all walls are the same, 
(b): Estimated critical exponents - a: -0.256, walls are wet, 
6: -0.262, walls are wet, c: -0.269, walls are hydrophilic. 

Thus, unlike earlier works, our results naturally account 
for the non-uniformity of the fluid and, due to the ab 
initio nature of our model which is based on statistical 
mechanics, can serve as "experimental" verification of the 
analysis in Ref. [I]. 

We solve the integral equations arising from the min- 
imization of |2j on a non-uniform collocation grid al- 
lowing us to capture the smooth decay of the fluid den- 
sity into the capillary bulk and its asymptotic behavior, 
which then determines the critical exponents. Note that 
we chose the interval of x < 100cr, for the grid to be suf- 
ficiently dense to resolve the transition of p (x, y) from 
capillary-liquid to capillary-vapor anywhere within that 
interval, even for temperatures far from T cc . The bound- 
ary condition of contact with the ID slit pore has been set 
up at the distance of x ~ 10 a. The numerical scheme 
has allowed us to follow the propagation of the liquid 
meniscus into the capillary bulk as A/i — > A^ c for a 
range of temperatures, wall potentials and capillary sizes 
(note that reducing H increases fluctuation effects, limit- 
ing the applicability of a mean-field treatment). Figure [3] 
summarizes our findings regarding the critical exponent 
for the divergence of T, which confirms well the analytic 
results of Parry et al. jj: T ~ (/j, — fic) 1 ^ 4 , and com- 
pletes our picture of the second-order CC transition in 
the infinite capped capillary. 

Finally, we report a multitude of non-concave branches 
of the isotherms which we find when the walls 




FIG. 4. Typical density profiles corresponding to the saddle 
points of the grand free energy in |TJ. 

are hydrophobic. Figure [4] shows some vivid examples of 
the density profiles corresponding to such branches (sad- 
dle points of f2[/o], eq. Q) in a finite square capillary 
(H = L), warranting an analogy with pattern formation 
phenomena. In a dynamic system an evolving density 
profile is likely to be attracted to such states, albeit tem- 
porarily, giving rise to a rather rich behavior [5J. See 
also Supplemental Movie 3. The fluid configuration re- 
flects the interplay of repulsive wall-fluid and attractive 
fluid-fluid interactions. A better understanding can be 
achieved, e.g., within a diffusive dynamic formulation of 
DF, [T3], or the one including hydrodynamic effects, |14) . 

We believe that our results will motivate further stud- 
ies on the role of confinement in fluids. The new phase 
transition described here, in particular, is an experimen- 
tally realizable scenario of the interplay between intra- 
fluid and fluid-substrate potentials. The dynamic investi- 
gation of confined systems also promises to provide rather 
interesting results related to contact line motion and the 
resulting pattern formation. 
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